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ABSTRACT 

The details of ignition of Type la supernovae remain fuzzy, despite the importance of this input 
for any large-scale model of the final explosion. Here, we begin a process of understanding the 
ignition of these hotspots by examining the burning of one zone of material, and then investigate 
the ignition of a detonation due to rapid heating at single point. 

We numerically measure the ignition delay time for onset of burning in mixtures of degenerate 
material and provide fitting formula for conditions of relevance in the Type la problem. Using the 
neon abundance as a proxy for the white dwarf progenitor's metallicity, we then find that ignition 
times can decrease by ~ 20% with addition of even 5% of neon by mass. When temperature 
fluctuations that successfully kindle a region are very rare, such a reduction in ignition time can 
increase the probability of ignition by orders of magnitude. If the neon comes largely at the 
expense of carbon, a similar increase in the ignition time can occur. 

We then consider the ignition of a detonation by an explosive energy input in one localized 
zone, e.g. a Sedov blast wave leading to a shock-ignited detonation. Building on previous work on 
curved detonations, we confirm that surprisingly large inputs of energy are required to successfully 
launch a detonation, leading to required matchheads of ps 4500 detonation thicknesses - tens of 
centimeters to hundreds of meters - which is orders of magnitude larger than naive considerations 
might suggest. This is a very difficult constraint to meet for some pictures of a deflagration-to- 
detonation transition, such as a Zel'dovich gradient mechanism ignition in the distributed burning 
regime. 

Subject headings: supernovae: general — white dwarfs - hydrodynamics — nuclear reactions, 
nucleosynthesis, abundances — methods: numerical 



1. INTRODUCTION 

The current favored model for Type la supernovae (SNIa) involves burning beginning as a subsonic 
deflagration near the central region of a Chandrasekhar-mass white dwarf. Progress has been made in 
recent years in understanding the middle stages of these events through multiscale reactive flow simulations 
where the initial burning is prescribed as an initial condition of one or more sizable bubbles already burning 
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material at time zero. However, the initial ignition process by which such bubbles begin burning - whether 
enormous 50 km bubbles (Plewa et al. 2004) or more physically motivated smaller igniting points (Rcincckc 
ct al. 2002; Hofiich and Stein 2002; Bravo and Garcfa-Senz 2003) remains poorly understood. Further, if 
later in the evolution there is a transition to a detonation (e.g., Gamezo et al. 2004), this ignition process, 
too, must be explained. Indeed, ignition physics will play a role — by determining the location, number, 
and sizes of the first burning points — in any currently viable SNIa model. However, until very recently (for 
example, Woosley et al. 2004; Garcfa-Senz and Bravo 2005) very little work has gone into examining the 
ignition physics of these events. Here we begin examining the ignition process by considering the simplest 
ignitions possible - that of a single zone - and the possibility of igniting a detonation from a Sedov blast 
wave launched at a single point. 

1.1. Ignition Delay Times 

Astrophysical combustion, like most combustion (for example, Williams 1985; Glassman 1996), is highly 
temperature-dependent; the 12 C + 12 C reaction, for example, scales as T 12 near 10 9 K. Rates for the 
exothermic reactions which define the burning process are generally exponential or near-exponential in 
temperature (e.g., Caughlan and Fowler 1988). Thus a region with a positive temperature perturbation can 
sit 'simmering' for a very long time, initially only very slowly consuming fuel and increasing its temperature 
as an exponential runaway occurs. This is especially true in the electron-degenerate environment of a white 
dwarf, where the small increases in temperature that occur for most of the evolution of the hotspot will have 
only extremely small hydrodynamic effects. 

If fuel depletion and hydrodynamical effects were ignored, the temperature of the spot would become 
infinite after a finite period of time. This time is called the ignition time, or ignition delay time, or sometimes 
induction time, Tj. After ignition starts, burning proceeds for some length of time t\>. 

For burning problems of interest, of course, fuel depletion is important, and no quantities become infinite; 
however, the idea of an ignition delay time still holds (see Fig. 1). If the energy release rate for most of the 
evolution of the burning is too small to have significant hydrodynamical effects, and if the timescale over 
which burning 'suddenly turns on' is much shorter than any other hydrodynamical or conductive timescales, 
then the burning of such a hotspot can be treated, as an excellent approximation, as a step function where 
all energy is released from t = Ti to t = r< + 75. In many problems, where Tb -C Tj, this can be further 
simplified to burning occurring only at t = n. Where such an approximation (often called 'high activation- 
energy asymptotics') holds, it greatly simplifies many problems of burning or ignition, reducing the region of 
burning in a flame to an infinitcsimally thin 'flamelet' (Matalon and Matkowsky 1982) surface, for instance, 
or the structure of a detonation to a 'square wave' (Erpenbeck 1963). Where this approximation does not 
hold - such as if slow /3-decay processes are important as bottlenecks for reactions to proceed (e.g., p-p 
burning or the CNO cycle) the simplification of burning happening only over n < t < Ti + t^ often remains 
useful. 

Even for the simple case of one zone, ignition delay times are relevant for investigation of ignition in 
SNIa because it sets a minimum time scale over which an initial local positive temperature perturbation (hot 
spot) can successfully ignite and launch a combustion wave; other timescales, such as turbulent disruption 
of the hotspot, or diffusive timescales, must be larger than this for ignition to successfully occur. 

If the burning occurs in an ideal gas, or in a material with some other simple equation of state, it is fairly 
easy to write down approximate ignition times for various burning laws. In a white dwarf, however, where 
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Fig. 1. — Temperature evolution for burning a zone at constant pressure with an initial state of Xi2 C = 1.0, 
T = 10 9 K, p = 5 x 10 8 g cm~ 3 . Because of the strong temperature dependence, a runaway takes place and 
most of the burning happens 'all at once'. 



the material is partially degenerate or relativistic and the equation of state is quite complicated (Timmes 
and Swesty 2000), no such closed- form expression can be written. In §2 we numerically follow the abundance 
and thermodynamic evolution of a zone of white dwarf material in order to measure the ignition times as a 
function of the initial temperature, density, and composition. We follow both constant-density and constant- 
pressure trajectories. The results are summarized by simple, moderately accurate, fitting formula. In §3 
we consider the ignition of a detonation through a localized energy release producing a Sedov blast wave, 
and estimate the amount of energy that must be released for the detonation to successfully ignite. In §4 
we consider our results in light of likely temperature fluctuation spectra during the simmering, convective 
phase. 



2. ONE ZONE IGNITION TIMES 
2.1. Calculations 

We performed a series of 1-zone calculations for the purposes of measuring ignition times in carbon- 
oxygen mixtures. For each of these two burning conditions - burning at constant volume and constant 
pressure - over 3500 initial conditions were examined, in a grid of initial densities, temperatures, and initial 
carbon fraction. The carbon (mass) fractions were in the range 0.4 < Xi2 C < 1.0, with the remainder being 
oxygen (Xwq = 1 — Xi2 C ), temperatures of 0.5 < Tg < 7, and densities 0.1 < p$ < 50, where Tg is the 
temperature in units of 10 9 K and p$ is the density in units of 10 8 g cm -3 . For burning, a 13-isotope a 
chain was used (Timmes 1999), and a Helmholtz free energy based stellar equation of state (Timmes and 
Swesty 2000) maintained the thermodynamic state. The temperature and abundance evolution equations 
were integrated together consistently, as described in Appendix A. Time evolutions were generated as in 
Fig. 1. 

To cover the wide range of burning times within each integration, the timestep was increased or decreased 
depending on the rate of abundance, energy release, and thermodynamic changes. The time interval was 
varied by up to a factor of two in each timestep, to try to keep the amount of energy released through 
burning per timestep change of fuel within the range 10~ 5 — 10~ 7 . Timesteps failing this criterion were 
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undone and re-taken with a smaller time interval. The final ignition time, when the simulation was stopped, 
was defined to be time when 90% of the carbon was consumed, although the time reported was found to be 
insensitive to endpoint chosen. The code used in this integration, as well as the resulting data, is available 
at http://www.cita.utoronto.ca/~ljdursi/ignition/. 

Over the initial conditions chosen, ignition times varied from 1CP 14 s to 10 +8 s. A representative contour 
plot showing the calculated ignition delay times are shown in Fig. 2. 
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Fig. 2. — Contour plot of ignition time as a function of initial density and temperature for a constant-pressure 
ignition of a mixture of half-carbon, half-oxygen by mass. 

Fitting formula for the ignition delay times under the two burning conditions were determined; the 
ignition time for constant-pressure ignition can be given approximately as 

n, cp (p, T, Xi2 C ) = 1.15 x 1(T 5 sec {X^p^ fc P (T) (1 + 1193/ cp (T)) (1) 
/ cp (T) = (T 9 - 0.214)- 7 ' 566 

and that for constant-volume ignition 

r 4 , cv (p, T, Xi2 C ) = 1.81 x 10" 5 sec (X^cPb)' 1 - 86 / CV (T) (1 + 1178/ CV (T)) (2) 
/cv(T) = (T 9 - 0.206)- 7 700 
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where T 9 is the initial temperature in units of 10 9 K, and p$ is the initial density in units of 10 8 g cm 3 . 

The fits are good to within a factor of five between 10 -9 sec and 1 sec, and to within a factor of 
10 between 10 -14 sec and 100 sec. This can be compared to other analytic expressions, for instance the 
constant-pressure formula from Woosley et al. (2004), 

T,=i5sec (A) 22 (») 3S (3) 

which, as shown in Fig. 4, is an excellent approximation over a somewhat more narrow range of conditions. 




Fig. 3. Fit results vs. calculated results for constant-pressure (left) and constant-volume (right) ignition 
delay times, for the full range of densities (0.1 < ps < 50), temperatures (0.5 < T 9 < 7), and carbon mass 
abundances (0.4 < Xi2 C < 1.0) considered. 




Fig. 4. — Woosley et al. (2004) ignition time results vs. calculated results for constant-pressure ignition delay 
times for a mixture half-carbon half-oxygen by mass (Xi2 C = 0.5), over a truncated range (1 < p$ < 50), 
(0.5 < T 9 < 1.5). 

Which of the two cases (constant volume or constant pressure) are appropriate will depend on comparing 
the ignition time with the relevant hydrodynamical timescale — in particular, the sound-crossing time of the 
region undergoing ignition. If the region is small enough that many sound crossing times occur during the 
ignition, then the constant-pressure value is appropriate; if ignition occurs in much less than a crossing time, 
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the constant-volume value is appropriate; intermediate timescales will result in intermediate values. In the 
case of turbulent ignition of a flame, presumably it is the constant-pressure time which will be most relevant. 
In any case, due to the degeneracy of the material in the density and temperature ranges considered here, 
the computed ignition times (or the fits) rarely differ between the two cases more than 50%. 



2.2. Effect of Metallicity 

Most of the initial metallicity of main-sequence star comes from the CNO and 56 Fe nuclei inherited from 
its ambient interstellar medium at birth. The slowest step in the hydrogen-burning CNO cycle is proton 
capture onto 14 N. This results in all the CNO catalysts piling up into 14 N when hydrogen burning on the 
main sequence is completed. During helium burning, all of the 14 N is converted into 22 Ne. 

As a proxy for investigating the effects of metallicity in our a-chain based reaction network, we consider 
the ignition of a Xi2 C = 0.5 constant-pressure ignition while increasing the fraction of 20 Ne (and thus 
decreasing the abundance of oxygen). We'll verify this surrogate by using 22 Ne in larger networks. 

The effects of increasing X2t> Ne from to 0.05, and then further to 0.1 and 0.2, is shown in Fig. 5. 
Addition of even fairly modest amounts of neon can significantly (20-30%) reduce the ignition times for 
much of the thermodynamic conditions evaluated here. The reduced ignition times are found to result from 
a larger energy release from burning over the entire integration range. 

To understand the increase in the energy deposited by burning, consider the a-chain reactions which 
dominate burning in this regime and are modeled by the 'aprox!3' network as: 



12 C + 12 C 20 Nc + a + 13.933 MeV 

12 C + a 16 + 7.162 MeV 

16 + a 20 Nc + 4.734 MeV (4) 

20 Ne + a 24 Mg + 9.312 MeV 

24 Mg + a 28 Si + 9.984 MeV. 

Given an initial mixture of 12 C, 16 0, and 20 Ne, it is the carbon burning reaction 12 C + 12 C which happens 
first. The resulting a-particle can capture onto carbon or any heavier element in the a chain. Unless that 
chain is already populated, then the (very exothermic) neon capture and all flows to still heavier nuclei are 
choked off until enough 20 Ne and 24 Mg are generated through burning. Adding even quite modest amounts 
of neon to the initial mixture allows more a-chain reactions to promptly occur during carbon burning. 

The plots in Fig. 6 show the abundance evolution of a-chain elements during constant-density burning 
at p$ — 10, Tg = 1, Xi2 C — 0.5, and an initial neon abundance of zero and 0.05. The inclusion of 5% neon 
by mass greatly speeds the production of heavier intermediate-mass elements, and thus the exothermicity of 
the burning, reducing the ignition time. 

To confirm that this effect is not artificially enhanced by using an a-chain reaction network, and to 
quantitatively verify the ignition times produced by the aproxl3 network used here, we compared the ignition 
times at two (p, T) points and varying neon abundance with those produced by reaction networks containing 
513 and 3304 isotopes. The results are shown in Table 1 We see that not only are the times computed 
with the smaller network quantitatively in good agreement with the larger networks (within 10% for the 
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Fig. 5. — The difference in ignition time for a constant-pressure ignition of Xi2 C = 0.5, Xie = 0.5 when 
some of the Oxygen is replaced by Neon-20. On the top left is shown the fractional difference in ignition 
time with the addition of X2o- Nc = 0.05 as a function of the ignition time with X2o Ne = 0.0. On the top 
right, bottom left, and bottom right are contour plots in p — T space of the percent difference in ignition 
time with Xm No = 0.05, 0.1, 0.2, respectively. The ignition time for the base case is shown in Fig. 2. 

lower-temperature case, and within 25% for the higher case), but the trends are similar. The same trend 
is apparent when 22 Ne, unavailable in the a-chain network, is used instead of 20 Ne. The trend is in fact 
stronger with 22 Ne because of additional flow paths that become available. 

While the above has demonstrated the active role neon plays in the ignition process, its effect is much 
smaller than that of carbon, which is the primary source of fuel in this burning process. As a result, in the 
more realistic case when an increase in neon comes at the expense of both carbon and oxygen, this effect 
is reduced (for small neon fractions) or completely reversed (for moderate neon fractions). This is shown, 
for instance, in Fig. 7. In this case, increased metallicity of the progenitor system has the opposite effect; it 
makes the ignition time significantly longer for hotspots in the resulting white dwarf. 
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Fig. 6. — Mass fraction evolution in constant-density burning, with an initial state of p% = 10, Tg = 1, 
Xi2 C — 0.5, with X20N e = 0.0 (left) and -X"2o Nc = 0.1 (right). The burning here was calculated with a 513- 
isotope network. Because of the removal of the a-chain bottleneck at neon on the right, burning proceeds 
faster and the generation of higher intermediate-mass elements is raised by orders of magnitude at early 
times. 



3. APPLICATION: IGNITION OF A DETONATION 

3.1. Detonation Structure 

Detonation waves are a supersonic mode of propagating combustion. A shock wave heats up material, 
which then ignites, releasing energy which further powers the shock. (See, for instance, textbooks such as 
Glassman 1996; Williams 1985). There are, broadly, four states in a detonation: the unshockcd material; 
the shocked material immediately behind the shock; an induction zone, where the heated material slowly 
begins burning and then the reaction zone, where the bulk of the exothermic burning takes place. 

Unsupported, self-sustaining detonations can be of the Chapman- Jouget (CJ) type, where at the end 
of the reaction zone the flow becomes sonic, or of the pathological type, where the sonic point occurs within 
reaction zone, decoupling the flow downstream of the sonic point from the shock. Pathological detonations 
can occur in material where there are endothermic reactions or other dissipative or cooling effects, and 
may have speeds slightly higher (typically by a few percent) than the CJ speed. Detonations within highly 
degenerate white-dwarf material (ps > 0.2) are of the pathological type (Khokhlov 1989), largely because 
some regions of the flow behind the detonation can have significant amounts of endothermic reactions, 
violating the assumptions of the CJ structure. Because in the case of a pathological detonation some 
fraction of the reactions powering the detonation are decoupled from the shock, calculating the speed of a 
pathological detonation requires detailed integration of the detonation structure, rather than simply using 
jump conditions. For cither kind of detonation, one can estimate where most burning occurs with shock 
speed D (which, even in non-CJ case, can be estimated with CJ speed) and n; k — Dri is the position 
behind the shock at which the induction zone ends and rapid burning takes place. 

In the case of a detonation into a very low-density, cold material, the material immediately behind the 
shock will still not burn significantly until a length of time equal to the ignition time passes, resulting in a 
'square wave' detonation. For conditions relevant to near the core of a white dwarf, however, the post-shock 
fluid will typically have temperatures on order Tg « 5 - that is, temperatures which are already near the 
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Fig. 7. — As in Fig. 5, the difference in ignition time for a constant-pressure ignition of Xi2 C = 0.5, 
Xi&q = 0.5 when (left) a mass fraction of 0.005 of each of the carbon and oxygen is replaced by Neon-20 and 
(right) when 0.05 of each is replaced by Neon. Note that in this case, for X2o Ne = 0.1, with the exception of 
a small region in (p,T) the ignition time is increased by the same magnitude that it is decreased in the case 
when all of the neon comes from oxygen, e.g. the Xi2 C = 0.5,Aie o = 0.4,X2o No = 0.1 case of Fig 5 in the 
bottom right panel. 

maximum temperature which will be obtained by burning. Even in these cases, this estimate of k provides 
a good measure of the thickness of the detonation structure behind the shock, as is shown in Fig 8 where it 
measures the position of maximum burning. 



3.2. Ignition of a Spherical Detonation 

One mechanism for ignition of a detonation by a hotspot is an initial rapid input of energy which leads 
to a Sedov blast wave; the material shocked by the outgoing spherical blast ignites, and as the outgoing 
wave slows, a steady outgoing detonation results. A steady detonation cannot form until the outgoing shock 
wave speed drops to the detonation velocity detonation, or else the energy released by reactions behind 
the blast wave will not be able to 'catch up' to the outgoing shock to drive it. On the other hand, if the 
shock drops significantly below the detonation speed before ignition takes place, not enough material will be 
burning per unit time to sustain a detonation wave. We denote the position of the shock when it reaches the 
detonation speed from above as Rd- Naively, the condition for successful detonation ignition would be that 
Rd ^ h, since a detonation structure with width of order k must be set up before the shock speed becomes 
too slow; a more sophisticated derivation of this criterion can be found in Zeldovich et al. (1956). However, 
experimentally this is known to be far too lenient a condition for terrestrial detonations and Rjy must be 
orders of magnitude larger than k (Desbordes 1986). This has also been empirically seen in the context of 
astrophysical detonations (e.g., Niemeyer and Woosley 1997). 

This has been explained by, for instance, He and Clavin (1994). Curvature has a significant nonlinear 
effect on the structure of a detonation — even more so than on the structure of a flame (e.g., Dursi et al. 2003) 
because the curvature directly effects the burning region rather than merely the preconditioning (diffusion) 
region. He and Clavin (1994), looking at a pseudo-steady calculation of a detonation with curvature, find 
that for a steady detonation to exist requires the curvature to be extremely small. The condition found 
by the authors requires Rd ~ 447 2 /(7 2 — l)/3zk, where 7 is the polytropic index of the ideal gas and (3z 
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Fig. 8. — Two examples of estimating the detonation thickness, U = v s Ti in a detonation. On the left, 
plotted energy release rate from nuclear reactions behind the shock of a leftward-traveling ZND detonation 
into a pure-carbon quiescent medium of p s = 1, T 9 = 0.05. The shocked state is ps = 2.97, Tg = 4.2, and 
the incoming fuel velocity behind the shock is 4.0 x 10 8 cm s^ 1 . For the shocked material, the predicted 
ignition time is « 3 x 10~ n sec. Even in this case, where the shocked temperature is so high that significant 
burning occurs immediately, and the 'square wave' detonation structure does not apply, the predicted = 
1.2 x 10~ 2 cm correctly matches the peak of the reaction zone. On the right, pressure (top) and energy release 
rate (bottom), plotted relative to their maximum values, behind a leftward-traveling slightly overdriven 
detonation into the same material as in the previous figure, calculated by the hydrodynamics code Flash 
(Fryxell et al. 2000; Calder et al. 2002). The line above the plotted quantities shows the predicted k calculated 
with the observed values in the shocked state. 

represents the temperature sensitivity of the burning law; for the reactions and temperatures of interest in 
astrophysical combustion, flz ~ 10 — 15 (Dursi et al. 2003). Using 7 = 4/3 to describe the highly degenerate 
material near the centre of the white dwarf, this would give Rd ~ 1000 — 1500^. 

While in many cases using a polytropic ideal gas equation of state to describe degenerate white dwarf 
material can be an excellent approximation, in combustion phenomenon where the burning rate is highly 
temperature dependent it is problematic, making the above result unsuitable. Further, the authors assume a 
CJ detonation, as opposed to the pathological detonation that occurs at high densities in a white dwarf. This 
makes a straightforward application of the results of He and Clavin (1994) to our problem of interest difficult. 
Given the complexity of the partially degenerate equation of state in the white dwarf, re-deriving the analytic 
results for this application would be difficult. However, the detailed effect of curvature on detonations in 
white dwarfs has already been studied in a different manner, by Sharpe (2001). In this approach, the velocity 
eigenvalue for the detonation - which depends sensitively on the detonation structure - is calculated using 
a shooting method to numerically integrate the structure of a pathological detonation to the sonic point. 
A range of possible detonation velocities is input, and then repeatedly bisected as the detonation structure 
with a fixed given curvature term is integrated assuming the current detonation velocity. The result is 
the velocity-curvature relation for a detonation into a given ambient material, and as a byproduct, the 
relationship between the thermodynamic structure (at least up to the sonic point) and the curvature for the 
curved detonation. 

We follow the method of Sharpe (2001), also described in Appendix B, and measure the maximum 
sustainable spherical curvature of a detonation using the same equation of state and burning network as 
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used in the ignition time study. An example of the detonation-speed versus curvature relation is given in 
Fig. 9. Beyond some maximum curvature K max , no steady-state unsupported detonation can exist; thus in 
the case of a spherical detonation, for a steady detonation to successfully ignite, Rd > t m ^. 
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Fig. 9. — Example of detonation speed vs. curvature, for a detonation into a quiescent medium of ps = 1, 
T 9 = 0.05, Xi2 C = 0.5, Xie = 0.5. 

We calculated the detonation speed versus curvature relation for 0.5 < p$ < 20 and 0.25 < Xi2 C < 
1.0. The unshocked material was set to a temperature of Tg = 0.05, although the results were seen to 
be insensitive to this value. Our results are shown in Fig. 10 and Table 2. The code used to perform 
the calculations is available at http://www.cita. utoronto.ca/^ljdursi/ignition. The estimated detonation 
thickness and a comparison with K max is given in Table 3. As compared with the He and Clavin (1994) 
results of R D > 1000 - 1500Z;, we find R D > 3000 - 6000^. 

We can approximately summarize our results for the detonation velocity and the maximum curvature 
of these curved detonations: 

( Y \ 0.0273 
^Y) Ps' 0678 (5) 

, x ^2.869 

Kmax = 3.128 x 10- 3 cm" 1 ( J p\- 272 (6) 
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where D is the planar detonation velocity, and K max is the maximum sustainable curvature. The fit for 
the maximum curvature is shown in Fig. 11. The velocity fits produce errors of less than 10% within the 
range examined; the fit for K max is better for highly degenerate material (within 25% for p 8 > 1) than for 
less-degenerate material (within a factor of 2.5 for the entire range considered.) 

One can estimate how much energy would have to be released at a point to produce a shock wave that 
slowed to the detonation speed D at a radius of Rd — «max by using the similarity relations of a Sedov blast 
wave propagating into a constant density material. The position of the shock is given by (Sedov 1959): 



where /3 depends on the equation of state. We performed numerical experiments for a Sedov blast wave 
propagating through degenerate material at p s = l,Tg = 0.05 with no burning and found /3 ~ 0.77. Given 
this relation, the shock velocity would be 



where e 32 is the point energy input in units of 10 32 erg, and ro is the radius in units of 1 cm. 

Our numerical experiments for calculating the (3 for the Sedov blast wave under these conditions were 
as follows. As described later in this section for the calculation of spherical shock-ignited detonations, 
we performed a series of one-dimensional spherically symmetric hydrodynamical simulations, in this case 
without any burning, for the hydrodynamic state described. We used very high resolution, with AMR 
(Ax k2x 10~ 4 cm, in a simulation box of 5 — 50cm for these simulations without burning) to ensure that the 
shock front was described adequately, and we placed the initial excess energy uniformly in the center-most 
eight zones (e.g., over 1.5 x 10~ 3 cm). We found that this combination of high resolution and finite size 
of energy input worked very well to produce the correct blast wave structure. We then tested an input of 
energy of 10 26 , 10 27 , 10 28 , and 10 29 erg, and for all input energies obtained a very clean Sedov-like scaling 
between position (or velocity) and time, with a scaling coefficient fit to be 0.77. As it turns out, these extra 
calculations were likely unnecessary; the Sedov blast wave simulations with burning also reproduce the same 
Sedov scaling with the same (3 until the speed of the blast wave slows to within a factor two of the planar 
detonation speed. 

The requirement that the detonation be successful, that is v s w v max at r > k^xi gives a minimum 
energy 



which for j3 = 0.77, p 8 = 1, Xi2 C = 1 gives e 32 = 0.18. 

We can test the applicability of this energy criterion to detonations in degenerate white dwarf material 
by running a reactive Sedov blast wave through a constant density material with an given input point 
energy and watch the success or failure of detonation ignition. We performed this calculation with the 
Flash code (Fryxell et al. 2000; Calder et al. 2002), and the results are shown in Fig. 12. The detonation 
was resolved with a finest resolution of Ax « 2 x 10~ 4 cm (e.g., 60 points within li w 1.2 x 10~ 2 cm) 
in a one-dimensional spherically symmetric domain of maximum radius 512cm. A very high resolution 
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was used to ensure the detonation structure was adequately resolved; improperly resolving detonation can 
greatly exaggerate curvature effects (Menikoff et al. 1996) as well as causing spurious ignition or incorrect 
propagation speeds (Fryxell et al. 1989). At time zero, the energy initial energy that causes the blast wave 
was deposited into a region of size 0.0015cm, or about eight zones; this is a large enough number of zones 
that the blast wave is well resolved even at early times, but a small enough region that the Sedov assumption 
of a point explosion remains valid, as demonstrated by the Sedov scaling behavior until very late times when 
burning significantly alters the shock solution. 

As shown in Fig. 12, successful detonations do indeed begin propagation at predicted the input energy 
but they are disrupted due to instabilities (e.g., Kriminski et al. 1998); this suggests that the energy re- 
quirement given here is a lower bound for successful detonation ignition. However, the results from these 
simulations need to be taken schematically rather than quantitatively. Because of the extremely large point 
energies required to successfully ignite a detonation, for much of the evolution shown in Fig. 12, the blast 
wave is moving superluminally, as no relativistic effects were included in the hydrodynamics! Part of this is 
an artifact of assuming that the energy input is at a very localized region. 

We can also consider the (more plausible) launch of a spherical shock and possible ignition of detonation 
from a region of finite size. Constraints on what that size must be will come from energetic considerations 
similar to that of the point explosion. If the required energy input comes predominantly from carbon burning 
(Ae w 5.6 x 10 17 erg g _1 ), then the volume of material which must be burned to provide such energy can 
be found by setting the minimum energy requirement from Eq. 10 to 4/37rr^pAe, or, plugging in numerical 
values, 2.09 x 10~ 7 (10 32 erg)/9 8 (Ai2 C /0.5)rk . This gives a required radius of burning region 



Note that in general this matchhead radius is actually larger than the minimum radius for a sustainable 
detonation. The two sizes are plotted as a function of density in Fig. 13. 

These results can be compared to empirical results from the astrophysics literature, for instance § 3.1.2 
of Nicmcycr and Woosley (1997), where the authors performed 100-zone spherically symmetric simulations 
to examine the scales of detonation. By imposing a linear temperature profile peaking at Tg = 3.2 over 
thirty zones of a given density and examining the smallest that region could still ignite a detonation. We 
compare our results, in particular k" 1 , to the size of their matchhcads. Some caution should be exercised 
in this comparison, as they are of different quantities — in our work we find the smallest possible radius of 
a steady-state detonation, whereas the quantity measured in Niemeyer and Woosley (1997) is the smallest 
size of a region of a particular imposed temperature profile which can launch a detonation. With that caveat 
in mind, a comparison of the results where they (nearly) overlap is given in Table 4. Because of the large 
amounts of energy imposed in the region by those authors, and because the temperature profile imposed 
implies that the detonation would normally 'begin' where the temperature drops to that of the ambient 
medium, the results are actually quite similar, with only the one point at Tg = 2 differing markedly. 
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4. 



DISCUSSIONS AND CONCLUSIONS 



Any currently feasible mechanism for the explosion of a Type la supernova involves propagating a 
burning front in a carbon-oxygen white dwarf. Correctly modeling the process which leads to ignition is 
crucial to understanding the resulting explosion; once a thermonuclear burning front begins propagating it 
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is very difficult to extinguish (Zingale et al. 2001; Niemeyer 1999) and differing locations, numbers and sizes 
of the ignition points can significantly change the large-scale character of the explosion (e.g., Garcfa-Senz 
and Bravo 2005; Plewa et al. 2004). 

In this paper, we have examined ignition at a single point by calculating and fitting ignition times 
over a range of thermodynamic conditions and abundances relevant for almost any mechanism of ignition 
for a Type la supernova. We have also considered the ignition of spherical detonations, determining sizes 
necessary for a detonation successfully ignite, and the energy required for a Sedov blast wave to lead to a 
shock-ignited detonation. 

Currently the favored model of SNIa ignition is for the central density and temperature to increase 
as a carbon-oxygen white dwarf grows by accretion from a main-sequence companion in a binary system. 
When the energy deposited by carbon burning exceeds neutrino losses, the white dwarf enters a simmering, 
convective stage that lasts 500-1000 years with convective velocities of order 20-100 km s _1 (e.<7., Hoflich and 
Stein 2002; Kuhlcn et al. 2005). Owing to the convection, the entropy of the core is almost constant. As 
the temperature continues to rise, thermonuclear flames are born at points with the highest temperatures 
at or near the center of the white dwarf. Hot, Ray lcigh- Taylor unstable bubbles begin to rise through the 
convective interior. 

In this turbulent convective picture of SNIa ignition, the first points to runaway will be rare events at 
the high-temperature tail of the turbulent distribution. Thus, it is important to accurately quantify ignition 
conditions. Our five-parameter fitting formulae (equations 2 and 3) reproduce our computed ignition delay 
times very well over 15 orders of magnitude, covering most thermodynamic conditions and compositions 
relevant for ignition in a carbon-oxygen white dwarf. In particular, we have found that increases in the 
white dwarf metallicity, modeled by increases in the 20 Ne abundance, can speed up the ignition delay time 
by ~ 30% if the neon comes largely from oxygen, or decrease the ignition delay time by a similar amount if 
it comes equally from carbon and oxygen. In the 'rare event' picture, if all else is equal, this suggests that 
the first ignition points may happen at significantly different conditions in metal-poor white dwarfs than 
in metal-rich progenitors - either the first ignition points could occur when the mean core temperature is 
somewhat lower, or more points could ignite. 

To quantify this effect, Fig. 14 shows the increase in the probability of a point igniting due to a 30% 
reduction in ignition time. Kuhlen et al. (2005) suggest that the distribution of temperature fluctuations 
during the simmering, convective phase is Gaussian, particularly at the positive-fluctuation end. Assuming 
this distribution and taking <5Trms as a free parameter that measures the strength of turbulence, one finds 
that if ignition is extremely rare, a 30% reduction in ignition time could produce a substantial increase in 
number in the (still rare) events. 

Another example where local ignition can change burning behavior is the propagation upwards through 
the star of a burning, Raylcigh- Taylor unstable bubble. Such a bubble will experience shear instabilities 
and shed sparks (Zingale et al. 2005). Whether these sparks ignite the unburned material into which they 
fall to produce more burning bubbles (or fizzle out) will change the burning volume, and thus the overall 
effectiveness of burning. However, the statistics of these sparks are not yet well characterized. 

In the second part of our work, we examine another aspect of point ignition - the possibility of igniting a 
detonation at a point. We have shown that very large amounts of energy are required to ignite a detonation, 
making the purely local ignition of a detonation extremely unlikely. In Fig. 13, we show the required 
minimum sizes for the successful ignition of a steady state detonation. The restriction on the curvature of 
a successful detonation restricts all models for ignition of a spherical detonation, for example placing limits 
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on the minimum size of a preconditioned region for successful ignition of a detonation by the Zel'dovich 
gradient mechanism (e.g., Khokhlov et al. 1997). 

In particular, one oft-mooted possibility for the transition to a detonation is for it to occur in the 
distributed burning regime near p% = 0.5. Multi-dimensional resolved calculations of burning in this regime 
(Bell et al. 2004) has found that, because of both the distributed burning itself and the resulting vigorous 
mixing of fuel and ash, the remaining pockets of unburned fuel can have Xi2 C as low as 0.1-0.2. Examining the 
case of Xi2 C = 0.125 and p$ — 0.5 we find that the restrictions on igniting a successful spherical detonation 
are especially stringent, requiring matchheads on order half a kilometer in size. On the face of it, this makes 
preconditioning of an ignition time gradient for the Zel'dovich mechanism seem fairly unlikely. 

In our consideration of one-zone ignitions, we have ignored the influence of hydrodynamics, which can 
transport energy into or away from a burning region. Further, in both the one-zone ignitions and our con- 
sideration of spherical detonations, we have ignored thermal transport which, in the electron-degenerate 
material of a white dwarf, can be rapid and efficient. It is certainly the case that in regions which do success- 
fully ignite, the ignition timescale will dominate that of hydrodynamics or thermal conductivity; however, 
without comparing this timescale to the competing timescales of thermal and hydrodynamic transport, we 
have no way of determining what the successfully igniting regions will be, and which regions will 'fizzle 
out'. Thus the results presented here can give, at best, one one piece of the story in determining ignition 
conditions. 

Comparisons between the ignition times (or shock-crossing time of the matchhead in the detonation 
case) to the hydrodynamic or thermal diffusion ignition times is relatively straightforward if a back of the 
envelope calculation is sufficient. Zingale et al. (2005) indicates that the buoyancy-driven turbulence relevant 
for ignition in SNIa follows a Kolmogorov turbulence, in which case the hydrodynamic time for destruction 
of a hot spot behaves as: 



where V CODV and L conv are the velocity and length scales of the large-scale convective motions that provide 
the stirring for the turbulence on the integral scales. 

Similarly, Woosley et al. (2004) gives approximate values for the electron-conduction-dominatcd thermal 
diffusivity relevant for the core of a white dwarf, for which one can find the diffusion time of a strong (several 
times the background temperature) Gaussian hotspot of size I: 



If one is in a regime when one of these timescales clearly 'wins', then this simple comparison of timescales 
is sufficient. However, in the case of the first ignition in the centre of a white dwarf, this will not be the 
case, as these timescales will necessarily be quite close. In this case, to understand the highly-nonlinear 
process of ignition, one must rely on computational experiments with all of the relevant physics included. In 
future work, we will consider the ignition of Gaussian hotspots in a quiescent medium with hydrodynamics 
and thermal diffusion, and then consider the addition of another piece of physics - the possibility of ignition 
under compression for the specific case of a pulsational delayed detonation. 
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Tabic 1. Constant-volume ignition times (sec) adding Neon with different reaction networks 



p s = lO,T 9 = l P8 = i,t 9 =3 

Xi2q Xi6q ^20n o ^22 Ne 13 isotopes a 513 isotopes 3304 isotopes 13 isotopes b 513 isotopes 



0.5 


0.5 


0.0 


0.0 


5.96 x 10" 


-2 


6.42 x 10~ 2 


6.42 x 10" 


-2 


2.16 x 10~ 8 


1.74 x 10~ 8 


0.5 


0.45 


0.05 


0.0 


5.77 x 10" 


-2 


6.49 x 10~ 2 






1.94 x 10~ 8 


1.58 x 10~ 8 


0.5 


0.4 


0.1 


0.0 


5.58 x 10" 


-2 


6.38 x to" 2 


6.38 x 10" 


-2 


1.79 x 10~ 8 


1.49 x 10~ 8 


0.5 


0.3 


0.2 


0.0 


5.18 x 10" 


-2 


6.09 x 10~ 2 






1.60 x 10~ 8 


1.37 x 10~ 8 


0.5 


0.45 


0.0 


0.05 






5.56 x 10~ 2 








9.61 x 10" 9 


0.5 


0.4 


0.0 


0.1 






5.54 x 10~ 2 


5.54 x 10" 


-2 




9.04 x 10" 9 


0.5 


0.3 


0.0 


0.2 






5.64 x 10~ 2 








8.92 x 10~ 9 



a Ignition time in seconds. The l aproxl3' network used primarily in this work gives ignition times within 10% of 
more complete networks for this (p, T), and with similar trends. 

b Ignition time in seconds. The l aproxl3' network used primarily in this work gives ignition times within 25% of 
more complete networks for this (p,T), and with similar trends. 
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Fig. 10. — Steady-state detonation velocity as a function of curvature for background densities of, from left 
to right and top to bottom, pg = (0-5, 1, 2, 5, 10, 20). In each plot, the velocities are shown for initial carbon 
abundances, top to bottom, of Xi2 C = (1, 0.75, 0.5, 0.25). So that they could be shown on the same horizontal 
scale, the curvatures have been scaled to the maximum sustainable curvature for each set of conditions, given 
by Table 2. Points represent calculated speeds, and solid lines are fits with parameters also given in Table 2. 
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Table 2. Measured and Fit Detonation Velocities as a function of curvature 









K b 
™max 




"max 






Fit Parameters 












(10 9 cm) 


cm -1 




10 9 cm 






b 


c 




d 




PS = 0.5 


Xl2 C = 


1.00 


1.207 


1.08 x 10" 


-2 


1.096 


-1.55 x 10" 


-2 


7.75 x lO" 1 


-1.26 x 10° 


1.51 x 10° 


Xl2 C = 


0.75 


1.181 


2.94 x 10" 


•3 


1.069 


-5.74 x 10- 


-4 


4.68 x lO" 1 


4.07 x 10- 


■ l 


1.26 x 10- 


-l 


Xl2 C = 


0.50 


1.152 


7.14 x 10" 


-4 


1.067 


-1.78 x 10- 


-3 


6.11 x lO" 1 


7.21 x 10- 


■ l 


-3.30 x 10" 


■ l 


Xl2Q = 


0.25 


1.120 


1.26 x 10" 


-4 


1.055 


-2.27 x 10" 


-3 


1.15 x 10° 


-1.39 x 10- 


■ l 


-7.54 x 10" 


■3 


Xl2Q = 


0.125 


1.110 


2.34 x 10" 


-5 


1.059 


-1.32 x 10- 


-2 


1.20 x 10° 


-1.65 x 10- 


■0 


1.37 x 10- 


•3 


PS = 1- 


Xl2 C = 


1.00 


1.218 


2.78 x 10" 


-2 


1.131 


-4.69 x 10" 


-3 


6.25 x 10- 1 


1.87 x 10- 


3 


3.77 x 10" 


-1 


Xl2 C = 


0.75 


1.195 


1.02 x 10" 


-2 


1.116 


-4.08 x 10" 


-3 


6.53 x lO" 1 


3.30 x 10- 


■1 


2.06 x 10" 


-2 


Xl2 C = 


0.50 


1.170 


3.15 x 10- 


-3 


1.102 


-5.41 x 10- 


-3 


8.30 x lO" 1 


2.64 x 10- 


■1 


-8.85 x 10" 


-2 


Xl2Q = 


0.25 


1.141 


5.58 x 10" 


-4 


1.085 


-2.99 x 10" 


-3 


1.17 x 10° 


-2.12 x 10- 


■1 


4.44 x 10- 


-2 


Xl2Q = 


0.125 


1.126 


8.90 x 10" 


-5 


1.092 


-9.94 x 10- 


-3 


9.55 x 10- 1 


1.35 x 10" 


■1 


-7.99 x 10- 


-2 


PS = 2. 


Xl2Q = 


1.00 


1.249 


6.67 x 10" 


-2 


1.176 


-5.52 x 10" 


-4 


7.93 x 10- 1 


8.76 x 10- 


2 


1.20 x 10" 


-1 


Xl2 C = 


0.75 


1.227 


2.72 x 10- 


-2 


1.169 


-2.72 x 10- 


-3 


6.75 x lO" 1 


4.44 x 10- 


■1 


-1.16 x 10- 


-1 


Xl2Q = 


0.50 


1.203 


8.90 x 10" 


-3 


1.152 


-1.07 x 10- 


-2 


9.50 x lO" 1 


1.07 x 10- 


■1 


-4.54 x 10- 


-2 


Xl2Q = 


0.25 


1.177 


1.56 x 10- 


-3 


1.133 


-8.23 x 10- 


-3 


1.10 x 10° 


-9.81 x 10- 


2 


2.26 x 10- 


-3 


Xl2Q = 


0.125 


1.163 


2.71 x 10" 


-4 


1.127 


-1.07 x 10- 


-3 


1.18 x 10° 


-2.05 x 10- 


■1 


3.07 x 10" 


-2 


Ps = 5. 


Xl2Q = 


1.00 


1.316 


1.96 x 10- 


-1 


1.257 


-5.81 x 10- 


-3 


8.76 x lO" 1 


1.20 x 10- 


1 


9.98 x 10" 


-3 


Xl2 C = 


0.75 


1.295 


8.49 x 10" 


-2 


1.241 


1.24 x 10- 


-3 


1.11 x 10° 


-2.25 x 10- 


■1 


1.11 x 10" 


-1 


Xl2 C = 


0.50 


1.273 


2.76 x 10" 


-2 


1.227 


-3.67 x 10- 


-3 


1.16 x 10° 


-2.27 x 10- 


■1 


6.94 x 10- 


-2 


X\2q = 


0.25 


1.248 


4.75 x 10- 


-3 


1.210 


-8.00 x 10- 


-3 


1.19 x 10° 


-2.25 x 10- 


■1 


4.18 x 10- 


-2 


Xl2Q = 


0.125 


1.235 


9.30 x 10" 


-4 


1.205 


-9.26 x 10" 


-3 


8.82 x 10- 1 


2.83 x 10- 


■1 


-1.56 x 10- 


-1 


PS = 10. 


Xl2 C = 


1.00 


1.386 


4.45 x 10- 


-1 


1.327 


-1.31 x 10- 


-4 


1.04 x 10° 


-1.53 x 10- 


1 


1.10 x 10- 


-1 


Xl2 C = 


0.75 


1.366 


1.94 x 10- 


-1 


1.316 


-2.65 x 10- 


-3 


1.08 x 10° 


-1.19 x 10- 


■1 


4.04 x 10- 


-2 


Xl2 C = 


0.50 


1.343 


6.27 x 10- 


-2 


1.298 


6.40 x 10- 


-3 


1.26 x 10° 


-3.40 x 10- 


■1 


7.62 x 10- 


-2 


Xl2 C = 


0.25 


1.319 


1.06 x 10- 


-2 


1.279 


6.21 x 10- 


-4 


1.40 x 10° 


-5.84 x 10- 


■1 


1.81 x 10- 


-1 


Xl2 C = 


0.125 


1.306 


2.24 x 10- 


-3 


1.269 


-5.93 x 10- 


-3 


1.35 x 10° 


-4.74 x 10- 


■1 


1.26 x 10- 


-1 


PS = 20. 


Xl2 C = 


1.00 


1.472 


1.03 x 10' 


J 


1.411 


1.03 x 10- 


-5 


1.14 x 10° 


-2.99 x 10- 


1 


1.59 x 10- 


-1 


Xl2 C = 


0.75 


1.451 


4.56 x 10" 


-1 


1.399 


-1.93 x 10- 


-3 


1.26 x 10° 


-4.25 x 10- 


■1 


1.68 x 10- 


-1 


Xl2 C = 


0.50 


1.429 


1.47 x 10" 


-1 


1.382 


2.49 x 10" 


-4 


1.41 x 10° 


-6.10 x 10- 


■1 


2.03 x 10" 


-1 


Xl2 C = 


0.25 


1.405 


2.41 x 10- 


-2 


1.363 


-4.72 x 10" 


-3 


1.43 x 10° 


-5.96 x 10" 


■1 


1.68 x 10" 


-1 



a The measured planar (k = 0) detonation speed. 

b The maximum curvature for which a steady-state self-sustained detonation can be found to persist. 
c The detonation speed at the maximum curvature K max 

d The detonation velocity as a function of k can be fit by v' = s] a + bn' + ck' 2 + dn' 3 , where v' = (v — V msix )/(D — 
fmax) and k' = (1 — re/ftmax). Note that for small curvature ft, one can define a 'Markstein length' (e.g., Dursi ct al. 
2003) of sorts by dv/dn = -((£> - Dmax)/(2Kmax))(& + 2c + 3d) /(a + b + c + d) 1 / 2 . 
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Table 3. Detonation Thicknesses compared with k. 



/ a K~ //• 



P8 = 0.5 



Xl2 C = 


1.00 


2.26 x 10- 


2 


4.09 x 10 3 


Xl2 C = 


0.75 


5.48 x 10- 


-2 


6.20 x 10 3 


Xl2 c = 


0.50 


1.78 x 10" 


-1 


7.85 x 10 3 


Xl2 C = 


0.25 


1.09 x 10° 


7.27 x 10 3 


P8 = 1- 


Xl2 c = 


1.00 


9.59 x 10" 


-3 


3.74 x 10 3 


Xl2 C = 


0.75 


2.31 x 10" 


-2 


4.23 x 10 3 


Xl2 C = 


0.50 


7.40 x 10" 


2 


4.29 x 10 3 


Xl2 c = 


0.25 


4.55 x 10" 


-1 


3.93 x 10 3 


P8 = 2. 


Xl2 C = 


1.00 


3.94 x 10" 


3 


3.81 x 10 3 


Xl2 C = 


0.75 


9.68 x 10" 


•3 


3.80 x 10 3 


Xl2 c = 


0.50 


3.16 x 10" 


-2 


3.55 x 10 3 


Xl2 G = 


0.25 


1.92 x 10" 


-1 


3.33 x 10 3 


ps = 5. 


Xl2 C = 


1.00 


1.17 x 10" 


-3 


4.35 x 10 3 


Xl2 C = 


0.75 


2.93 x 10" 


3 


4.02 x 10 3 


Xl2 C = 


0.50 


9.58 x 10" 


3 


3.78 x 10 3 


Xl2 C = 


0.25 


6.02 x 10" 


-2 


3.50 x 10 3 


P8 = 10. 


Xl2 C = 


1.00 


4.51 x 10" 


■4 


4.99 x 10 3 


Xl2 C = 


0.75 


1.12 x 10" 


-3 


4.58 x 10 3 


Xl2 C = 


0.50 


3.89 x 10" 


-3 


4.10 x 10 3 


Xl2 C = 


0.25 


2.47 x 10" 


-2 


3.82 x 10 3 


P8 = 20. 


Xl2 C = 


1.00 


1.71 x 10" 


-4 


5.67 x 10 3 


Xl2 C = 


0.75 


4.48 x 10" 


-4 


4.89 x 10 3 


Xl2 C = 


0.50 


1.55 x 10" 


3 


4.40 x 10 3 


Xl2 c = 


0.25 


1.01 x 10" 


-2 


4.09 x 10 3 



a Detonation thickness, in cm, estimated 
by ij = v s Ti, where v s is the velocity of 
the shocked material immediately behind the 
leading shock and n is the ignition time at 
the same position. 

b «max is taken from Table 2. 
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Fig. 11. — Fit of maximum curvature for a sustainable steady-state detonation; the fit is within 25% for 
ps > 1, and within a factor of 2.5 for the entire range. 
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0.1 1.0 10.0 100.0 1000.0 

Shock Position (cm) 

Fig. 12. — Shock velocity vs. shock position for a Sedov blast wave with burning, propagating through 
constant density material. The velocity is measured in units of the steady planar detonation velocity, 
D = 1.218 x 10 9 cm/s. Input energies are, left to right, 10 27 , 10 28 , 10 29 , 10 30 , 10 31 , 10 32 ergs. Note that when 
the shock velocity approaches the detonation velocity from above, burning begins to effect the shock speed - 
with transients robustly occurring when the shock speed falls to about twice the planar detonation velocity 
- but the shock fails to sustain a detonation for the case of input energies of less than 10 31 ergs. For the 
two largest energies, a detonation does successfully propagate for many detonation thicknesses - which, for 
comparison purposes, is about 10~ 2 cm - but is disrupted by instabilities. The inset shows a closeup of the 
highest-energy 'detonation' becoming unstable, although it propagates w 10 3 ^. 
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Fig. 13. — Minimum scales for detonation. Plotted is k~^ x from Table 2, with dashed lines and points, and n, 
from Eq. 11 with solid lines. An expending shock must slow to approximately the steady detonation velocity 
at a radius no less k^Lo an d r b is the radius in which is contained enough potential nuclear energy from 
carbon burning to produce the energy for a 'Sedov' explosion which would would naturally slow down to D at 
a radius greater than K~a X . ^ n both cases, lines are plotted for, top to bottom, Xi2 C = 1/8, 1/4, 1/2, 3/4, 1. 
For the cases considered here, r& > k~^x- 
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Table 4. Minimum scales for detonation ignition calculated here and in Niemeyer & Woosley (1997) 



P8 Xi2 C * Xi 60 a NW97 b thiswork c 



1 


1 





40 cm 


36 cm 


1 


i 

2 


1 

? 

2 


2 m 


3.17 m 


20 


1 

2 


70 cm 


6.8 cm 


0.3 


1 





1 m 


30.2 cm d 


0.3 


1 

2 


1 

2 


50 m 


14.8 m d 



a Mass fraction of Carbon-12 or Oxygen-16. 

b Size of region of imposed linear tempera- 
ture profile, with peak temperature at Tg = 
3.2, needed to provide a detonation in a 
100-zonc spherically symmetric simulation in 
Niemeyer and Woosley (1997). 

c ^max as measured in this section. 

d ^max extrapolated from measurements us- 
ing the fitting formula given in Eq. 6. 
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Fig. 14. — Plots showing enhancement of probability of ignition for a 30% reduction in ignition time for a 
fiducial case (Tg = 2, p$ = 10, Xi2 C = 0.5, and requiring ignition within 0.1ns). On the left is shown the 
base-10 logarithm of the probability of ignition as a function of the RMS temperature fluctuation (5Trms/T) 
with and without the 30% reduction in ignition time - e.g., log 10 P(ri iCp < 10~ 10 sec) and log 10 P(0.7ri. cp < 
10~ 10 sec) - for the constant pressure ignition time fit given in Eq. 2 and assuming Gaussian distribution 
of temperature fluctuations. On the right is shown the base-10 logarithm of the probability enhancement 
of ignition with the 30% reduction in ignition time, e.g. log 10 [-P(0.7Ti. C p < 10 -10 sec)/P(Ti. cp < 10~ 10 sec)]. 
When ignition is very rare (such as <5Trms < 0.1T for this case), the increase in probability of ignition with 
a 30% ignition temperature decrease can be substantial (2-7 orders of magnitude). 
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A. Constant-Pressure and Constant- Volume Burning 

In this appendix, we describe the coupling of temperature evolution and burning used in the reaction 
networks used in this work. 

For a Helmholtz free energy based thermodynamic system, the temperature T, density p, and molar 
composition vector Y are the primary variables. (Yi are related to the mass fractions used elsewhere 
in this paper through = AiY i: where A t is the atomic mass number of the i th element.) We'll want to 
expand thermodynamic quantities in terms of these variables, for example, 

dP dP dP 

dp -df dT + ^ + ^W/ Y - (A1) 

All reasonable stellar equations of state will return the dP/dT and dP/dp partial derivatives (for a 
survey, see Timmes and Arnett 1999). The last term in equation (Al) is a bit tricker. Nearly all stellar 
equations of state use the mean atomic weight A and the mean charge Z to characterize the composition. 
Thus, the last term in equation (Al) can be expanded as 

? m 1 = ? dim 1 + ? dzM 1 • (A2) 

ii % 

All complete stellar equations of state will return the dP/dA and dP/dZ terms, so we will not be concerned 
with those terms (e.g., Timmes and Swesty 2000). In terms of the molar composition the mean atomic 
weight A and the mean charge Z are 

A =ETn * = W- <A3) 



The partial derivatives of A and Z are then 



<9A / 1 



(A4) 



dZ / 1 \V„„ Z, 



J2 Z * Y * + 



= -A 2 Y / Z l Y t + A 2 Z t 



A^Z.Y^j 



= A\Z t 

= A(Zi-Z) (A5) 
surprisingly simple. Operating with d/dt on equation (Al) yields a pressure evolution equation 
dP dPdT dPdp 9Pv- dAdYi_ dPy- dZ^dYi_ 

The last two terms of equation (A6) are easy to compute since the dYi/dt are the right-hand sides of the 
ordinary differential equations that comprise a nuclear reaction network equation. 
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The First Law of Thermodynamics 



dE + Pd 



G 



operated on with d/dt in the presence of an energy source becomes 

dE P dp _ ^ 
~dt ~f?~dl ~ 

Expanding the specific internal energy dE/dt term 

dE dT dE dp dE ^ OA dY, dE ^ dZ dYi P dp _ ■ 
^^ + ^dt + M^W l lU + ^^W l ~dT~'l?dt^ 



(A7) 



(A8) 



(A9) 



and collecting like terms becomes 



dEdT 
df~dt 



dE 
~dp~ 



dp dA dY, dE ^ dZ dYi _ ■ 



dt dA ^ dY, dt dZ 



dYi dt 



(A10) 



For a constant density evolution, dp = 0, and equation (A10) may be written as 



dT 
~dt 



dE 
dT 



■ dE dA dY, dE dZ dYi 
" " ^ dY,~dt 



dA 



dZ 



dY dt 



(All) 



Equation (All) obeys the First Law of Thermodynamics while being consistent with a general stellar equation 
of state and the constraint of constant density (no PdV work). To self-consistently couple the thermodynamic 
and abundance evolutions we implicitly solve equation (All) with the constraint of a constant density and 
the reaction network equations 



d ± = o 
dt 

dY C 



p 2 N\Rj^YjY k Yi 



(A12) 
(A13) 



as a single, fully coupled system of ordinary differential equations. We use the semi-implicit Bader-Denhard 
algorithm with the MA28 sparse matrix package as described in Timmes (1999). to solve equations (All) 
and (A12) simultaneously. In equation (A12) the (stoichiometric) coefficients Ci can be positive or nega- 
tive numbers that specify how many particles of species i are created or destroyed. The factorials in the 
denominator avoid double counting. The Ri are the reaction rate for the species involved, with the first 
sum representing reactions involving a single nucleus, the second sum representing binary reactions, and the 
third sum representing three-particle processes. 



For a constant pressure evolution, dP = 0, equation (A6) becomes 



dP dT dP dp 
dT~dt + ~dp~~dt 



dP ^ dA dY 

E 



dA 



dP v-^ dZ dY, 
E 



dYi dt dZ ^ dYi dt 



(A14) 



Solving for dp/dt yields 



dp 
~dl 



dP 

dp 



dPdT aP^_9Arfy, dP_y,dZdYi 



dT dt dA 



dY, dt 



(A15) 
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Substituting equation (A15) into equation (A10) 



dE_dT_ 
dfddi 



' dE _ P_ ' 




dp p 2 




dP 




L dp . 





dPdT 
df~dt 



dP 



+ 



dE 
M 



E 



dA dY t 
B?i dt 



dE 



E 



dk dYj 
Wi~dT 

dZ dY t 
dYi dt 



Eu^ aii dP \ -. 
fiY~M~ + ~d~Z 



dZ dY t 
dYi dt 



(A16) 



= S 



and collecting like terms yields 



dE 
df 

dE 

~dA 

dE 



dE 


P 


dp 


I' 2 


dP 
dp 




dE 


P 


dp 


1? 


dP 




dp 




dE 


P 


dp 




dP 




dp 






dT 
dt 



+ 



(A17) 



x ^ dA dYj 
^dYi~df + 

s ^dZdY L= . 

^ dY t dt 



and solving for dT/dt yields 



dT 
~dt 




(A18) 



dA dY t 
dY,~dT 




^ dZ dY t 
^ dY, dt 



Equation (A18) obeys the First Law of Thermodynamics (including PdV work) while being consistent with 
a general stellar equation of state and the constraint of constant pressure. To self-consistently couple the 
thermodynamic and abundance evolutions we implicitly solve equation (A18) along with the constraint of a 
constant pressure and the reaction network equations 



9 
dt 



f{P Pqotis 



0) 



a 



E^y^^ + E 



Ci 



CAC k \d\ 



p 2 N\R^ k YjY k Yi 



(A19) 
(A20) 



as a single, fully coupled system of differential-algebraic equations. We use the same integration method and 
linear algebra solver listed above to solve equations (A18) and (A19) simultaneously. Given the constant 
upstream pressure, temperature, and composition, a root-find is performed to find the density p in equation 
A19 such that the pressure returned from an equation of state is equal to the constant upstream pressure. 
Note the similarities between equation (All) for constant density and equation (A18) for constant pressure. 
Despite its complex appearance, equation (A18) can be concisely coded. The method described here is 
implemented in the burning networks found at http://www.cococubed.com/code_pages/codes.shtml. 

A complementary approach for isobaric combustion evolution was recently formulated by Cabezon 
et al. (2004). Although their approach is tied to the first-order Euler integration method, they demonstrated 
satisfactory accuracy and consistency of the thermodynamic and abundance evolutions. A direct comparison 
of the two methods should be examined in future investigations. 
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B. Calculation of Steady-State Detonation Structures 

Here we describe our method, which is almost completely as described in Sharpe (2001), for calculating 
the steady state velocities for unsupported detonations under curvature. 

The starting point is to consider the quasi-steady equations of hydrodynamics in shock attached frame, 
assuming curvature is constant through structure. This gives us Sharpe's eqns (6)-(9): 



dpu n 
dn 



+ up (u n + D n ) = 



(Bl) 



du n 1 dp 

u n — h -— = 

dn p dn 



(B2) 



de 
dn 



p_dp_ 

p 2 dn 



(B3) 



dX t 
dn 



Ri 



(B4) 



where n is the coordinate normal to the shock, u n and D n are the fluid velocities relative to the shock 
and D n is the normal detonation velocity Xi is the mass abundance of the i th species, Ri is net production 
rate of the i th species, p is the gasdynamic pressure, and e is total specific energy of the fluid. 

Using arguments similar to those described in the previous appendix, one can construct ODEs for the 
post-shock density and temperature structure in the detonation: 



dp 
dn 



a) - ul 



dT 
dn 



dp 
df 



■> / dp\ 1 dp ^ / dp \ dXi \ 



(B5) 



(B6) 



where o/, the frozen sound speed, is given by 



2 ( d P 



T,Y 



8e\ 

dp) TY 



dp 
df 



de_ 



(B7) 



and <f>, the thermicity, is given by 

13 

9X lJ PiTiYjt jjti 




de 

dx, 



p,T,Y jt j^i 



dX l 



(B8) 



Note that at the sonic point - as u n a/ - the density becomes unbounded unless at the same time, 
<f> — > 0. This pair of conditions defines the generalized sonic point, and defines a required boundary value at 
that singular point. 
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Normally it would be best to integrate from the singular point to the shock that defines the detonation 
front, but we do not have enough initial conditions to do this. Instead we use a shooting method; given 
an input range of detonation velocities, we pick the middle value and integrate these equations using an 
IVODE solver until we hit the singular point. If the initial guess was close enough to the correct detonation 
velocity, then if we hit <f> = before the sonicity becomes zero, our guessed velocity was too high; if we 
hit u n = df first, the guessed velocity was too low. We then continue bisecting the range and choosing 
the new midpoint a fixed number of times until we have the velocity within a satisfactorily small range. 
(The velocities were calculated to ranges much smaller — 2-20cm s _1 — than suggested by the number of 
significant figures quoted in our results here). 

The only difference in how we proceeded compared to Sharpe (2001) concerned how the partial deriva- 
tives were taken with respect to Xi. Sharpe used a 13-isotope burning network, as do we, but did not evolve 
the 13th isotope, instead explicitly using the definition of mass fraction, 



Thus (presumably) differentiating by was performed computing finite differences by varying JQ, and 
implicitly changing A 13 by the same amount in the opposite direction. However, there are arbitrarily many 
other possibilities that could be used to 'redistribute' the change made to Xj. We found it more accurate 
to consider the changes to the other mass fractions to be proportional to the amount the abundances were 
varying along the detonation structure at that point; these values were available from the burning routines. 

The code used in computing these structures and velocities is available at http://www.cita.utoronto.ca/~ljdursi/ignition. 
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